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We study the non-Markovian dynamics of a qubit made up of a two-level atom interacting with 
an electromagnetic field (EMF) initially at finite temperature. Unlike most earlier studies where the 
bath is assumed to be fixed, we study the coherent evolution of the combined qubit- EMF system, thus 
allowing for the back-action from the bath on the qubit and the qubit on the bath in a self-consistent 
manner. In this way we can see the development of quantum correlations and entanglement between 
the system and its environment, and how that affects the decoherence and relaxation of the system. 
, We find non-exponential decay for both the diagonal and non-diagonal matrix elements of the qubit 's 

reduced density matrix in the pointer basis. From the diagonal elements we see the qubit relaxes to 
thermal equilibrium with the bath. From the non-diagonal elements, we see the decoherence rate 
beginning at the usually predicted thermal rate, but changing to the zero temperature decoherence 
rate as the qubit and bath become entangled. These two rates are comparable, as was shown before 
in the zero temperature case [C. Anastopoulos and B. L. Hu, Phys. Rev. A 62 (2000) 033821]. On 
the entanglement of a qubit with the EMF under this type of resonant coupling we calculated, for 
the qubit reduced density matrix, the fidelity and the von Neumann entropy, which is a measure 
£f>^ ' of the purity of the density matrix. The present more accurate non-Markovian calculations predict 

lower loss of fidelity and purity as compared with the Markovian results. Generally speaking, with 
the inclusion of quantum correlations between the qubit and its environment, the non-Markovian 
processes tend to slow down the drive of the system to equilibrium, prolonging the decoherence and 
^ ' better preserving the fidelity and purity of the system. 
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O ■ !• INTRODUCTION 

Interest in quantum entanglement has grown in recent years motivated mainly by the attempt to understand and 
Q ' realize quantum information processing. An important aspect of quantum entanglement, which is unavoidable in 
quantum information processing, is the entanglement of a system with its environment. This issue is especially im- 
portant to the feasibility of quantum computation, as the error threshold for fault tolerant error correction (without 
which quantum computation is impossible) , is highly sensitive to the degree to which the environment becomes corre- 
lated with the qubits pj . In many realistic quantum computing schemes the environment includes an electromagnetic 
qh] field (EMF) . Studies of entanglement between multi-partite systems exist 13, 0. [J 0. @ . but few have attempted to 
predict the effects on a qubit from its entanglement with the EMF 0, S 0> D3 • We have addressed the relaxation and 
, ' decoherence issues in various contexts, from a two-level atom in a (zero-temperature) EMF [TlT ] to moving atoms in 
a cavity (HEi. Here we continue this study for these two issues and the issue of entanglement for a two-level atom 
!-h ' in a finite temperature EMF. 

In this paper we aim at addressing the issue of system-environment entanglement by carefully analyzing its effect 
on reduced system dynamics. Specifically, we study the reduced dynamics of a single qubit interacting with an 
initially thermal bath in the multimode Jaynes-Cummings model. The model is a well studied example of open 
system dynamics, however . pr ior anal yses have focused on Markovian dynamics by assuming no disturbance of the 
bath modes by the qubit |Tj, [la, lid 1171 Il8l [i9| . Although such a simplifying assumption does include a level of 
back-action, it neglects entanglement that forms between the qubit and bath modes during the evolution. Earlier 
predictions based on simplified Markovian approximations are thus unsuitable for studies of such entanglement, and 
should be scrutinized carefully before being applied to quantum computation (see e.g. |20| for a discussion of this issue 
bearing on error correction) . More realistic physical conditions are better served by non-Markovian treatme nts, which 
have begun to appear. However, many of them use approximations with limited short-time validity [2U l22l 1231 l24j or 
unclear physical meaning such as in non-Markovian stochastic Schroedinger equations [25L |2(J . 
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Our approach includes bath dynamics as well as qubit dynamics, and the quantum correlations between them. 
Analytic expressions are derived for the qubit reduced density matrix elements at low temperature, from which the 
fidelity (defined as the overlap between the initial and evolved qubit state) and von Neumann entropy are computed. 
Inspection of the matrix elements themselves show slowed (sub-exponential) decoherence for the off-diagonal elements 
and slowed relaxation to thermal asymptotes for the diagonal matrix elements, as compared to Markovian predictions. 
Analysis of the fidelity and von Neumann entropy similarly show slowed loss of fidelity and purity in the case of non- 
Mar kovian dynamics over Markovian dynamics. The overall picture which develops is of increased preservation of 
coherence in non-Markovian dynamics. 

Following successive degradation from an exact solution to a Markovian description of system-bath interactions, 
three distinct approximations are usually invoked. They are 1) the 2 nd order Born approximation, 2) the I s * Markov 
approximation 19], and 3) the assumption of a bath which is unaffected by its interaction with the system. 1) The 
2 nd order Born approximation is an approximation in the strength of the coupling constant, and applying it neglects 
terms of higher than 2 nd order in the coupling [2]j. 2) The I s * Markov approximation is an approximation in the 
back-action correlation time. It is a Markov approximation because it assumes that the back-action of the system 
onto itself through the bath at time t will depend only on the state of the combined system-bath at time t, and not 
on the past history. It is called the I s * Markov approximation because although it depends only on the state of the 
system at time t, it depends on the state of the bath as well as the system (through the bath averages), thus including 
the bath dynamics 19]. 3) The last of the above three approximations is the assumption of a bath state which is 
fixed for all time. That assumption expressly excludes any dynamical evolution of the bath. 

In the usual derivation of the Markovian master equation from the Schrodinger equation for the system-bath 
density operator, all three approximations are made. In contrast, the Heisenberg-Langevin approaches make only 
the I s * Markov approximation. However, for spin-boson models such approaches have focused on the reduced qubit 
dynamics in strictly vacuum EMF, although in the presence of a classical source (e.g. resonance florescence). The 
resulting equations for the qubit degrees of freedom are called the Bloch-Langevin equations ■ The Schrodinger- 
master equation can be derived from the Heisenber g-L angevin equation after a perturbative expansion which imposes 
the first and third approximations from the above |28| . For a comparison of these approximations see Ref. [2?j . Our 
path integral approach to the reduced system dynamics uses only the first of the above three approximations by 
allowing the combined system-bath to evolve coherently throughout the interaction period. Only at the end of all 
coherent evolution will the bath variables be traced out to yield the reduced system evolution. 

The approach we take is straightforward, although the actual implementation includes some nonstandard techniques 
involving Grassmann path integrals. In Section II the Hamiltonian and other important aspects of the model, including 
the coherent state represenation, are reviewed. The transition amplitude is derived in Section IIIA, utilizing the 
coherent state representation for the bosonic degrees of freedom and Grassmann states for the qubit degrees of 
freedom, following 11]. Doing so will involve a recursive computation which exploits the semigroup property of 
the transition matrix (similar to a technique used in Ref. |30|). After evaluating the transition amplitudes in an 
intermediate form, we combine the forward and backward versions by tracing over the final bosonic coherent states to 
construct the reduced propagator in Section IIIB. An initial thermal state for the oscillator bath is then introduced 
and the reduced dynamics of the qubit are calculated in Section IIIC. Section IV gives discussion and further analysis 
of the results. 



II. MODEL 
A. Hamiltonian 

The model used for atom-field interactions is the standard multi-mode Jaynes-Cummings model of a two-level 
system interacting with a harmonic oscillator bath. The total Hamiltonian under the dipole, rotating wave and 
two- level approximations is given by (e.g. Appendix A of [TTj]) 

H = hu S z + hY, [ukb{b k + (\ k S + b k + hS-bfy] (1) 
k 

where bl,bk are the creation and annihilation operators for the k th bath mode with frequency uik, and hui is the 
energy separation between the two levels. The operators S z , S+, and S— are the qubit operators for z- projection, 
spin-up, and spin-down, respectively. The couplings, and have absorbed a dependence on the spectral density 
of the bath [13 . 
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B. Coherent States 



Coherent states are denned as any set of states generated by the exponentiated operation of a creation operator 
and a suitable label on a chosen fiducial state [^J |32| : 

\z k ) = cxp(z fc ^)| 0fc ) (2) 
\rj) = cxp(^+)|0) (3) 

In the case of the bosonic coherent states, defined in Eq. J2Jl, the label z k is a complex number, and in the case of the 
Grassmann coherent states, defined in Eq. the label r\ is an anti-commuting number. The chosen fiducial states 
are the EMF vacuum and the lower state of the two-level system, respectively. A state of the combined atom-field 
system can be expanded in a direct product coherent state basis, 

\{z k },ri) = \{z k })®\ri), (4) 

in which the bosonic coherent state basis, \z k ), is used to represent the EMF, and a Grassmann coherent state basis, 
|ry), is used to represent the two-level internal degrees of freedom of the atom. 

Grassmann and bosonic coherent states share well known properties of general coherent states, such as being 
non-orthogonal and eigenstates of the annihilator, 

(z k \z' k ) = exp(z k z' k ) (f?|V) = exp(fjr)') (5) 
b k \z k ) = z k \z k ) S-\rf} = r)\rj), (6) 

where the overbar denotes conjugation. Despite their non-orthogonality, both types of coherent states are (over- 
complete sets and have a resolution of unity, 



d/j,{z k )\z k )(z k \ = / d/z(7?)|77)(j7| (7) 



with the measures 



dfi(zk) = exp(-z k z k )dz k dz k (8) 
d/z(?j) = exp(— fjr])dfjdr]. (9) 

That these measures are exponential functions makes the coherent states a particularly suitable representation for 
transition amplitudes written as path integrals. For convenience the short hand notation 

]Jdfi{z k ) = dn{{z k }) (10) 

k 

is defined for the product of the measure of different mode coherent states. 

In the bosonic and Grassmann coherent states, the Q-representation of the Hamiltonian, Eq. QJ, is 

(fj, {zk}\H\r}', {z' k }) = huj fjr)' + [ u kZhZ k + (Xkf}z k + \ k z k -q')] , 

k 

in which the replacement S z — > S+S-, correct up to an additive constant, was made. The Q-representation Hamilto- 
nian will participate prominently in the path integrals of the next section. 



III. APPROACH 



A. Transition Amplitude 



Here we construct and evaluate the transition amplitude K(tf,ti) of coherent states from an initial time, tj = 0, to 
coherent states at a final time, tf — t, 

K(t,0) = (fjf,{z fk y t t\U(t,0)\rH,{z*hO)- ( n ) 
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with U(t, 0) being the time evolution operator, 

U(t,0) = e -iti Hds . (12) 

Following the path integral methodology, we partition the interval [0, t] into a large number (TV) of time steps, such 
that t — Ne. The path integral is then calculated as a discrete functional. Doing so, the n-step transition amplitude 
can be written in a general form, 



K(ne, 0) = exp <^ rj n il> n + ^ z nk f nk + ^ r] n g nk + ^ 

^ l, U I. J 



(13) 



k k k 

By applying the semigroup property of the transition amplitude, 



K({n + l)e, 0) = / dfi{ Vn ) J dfi({z k })K((n + l)e, ne)K(ne, 0) (14) 

Unite difference relations can be found for the coefficients in the action. Setting h = 1, and absorbing factors of 2n, 
they are, 

Ipn = (1 - «W e)V>n-l + Y^ki^nM^n-hk IpO = Vi Qcx 

4>n,k = («A n ,fee)^„_i + (1 - iuj k e)(/)n-i,k </>o,fc = 

9n,k = iw t)g n -\,k + (i^n,k€)fn-l,k 90,k = , . 

/n,fc = (iA n> fee) + (1 - iUke)fn-l,k fo,k = Zi,k- 

The coupling constants in Eqs. (|15I16|) have time indices because they are separated by complete sets of states 
at different time steps when the Hamiltonian is partitioned, thus they are separate sets of Grassmann pairs. The 
transition amplitude at time t(= Ne) can be written: 

K{t, 0) = exp^/VCt) + ^2 E fkfk{t) + ^2 Vf9k(t) + ZfMt)}. (17) 

Since this equation is a function of Grassmann variables it is to be treated as a formal expression that has meaning 
only in its polynomial expansion. In that polynomial expansion many terms will be truncated due to the nilpotency 
of the Grassmann variables. Expanding out Eq. Ijl7|l and defining the functionals 



F[{m e }](t) 


k 


(18) 


GiUm^t) 


= 9i(t)Y[(fk(t)) m * 

k 


(19) 


*'[{«*}](*) 


= m\{Uk{t)) mk 

k 


(20) 




= Mt)\{Uk{t)) mk 

k 


(21) 




= 9 P (t)Mt)Y[(fk(t)) mk 

k 


(22) 




= 9 P {t)m\[{fk{t)T k 


(23) 



gives the following expanded expression for the transition amplitude (with time dependence left implied for notational 
clarity) 



K(t,0)= £ 

{m £ } 



F[{m c }]+fj f ^ f [{m c }]+J2vfGi[{m c }] 

,24) 

E% $ p[^}] + Y l fjfZf P ®i p [{m s }} 
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The variable is the number of photons in the ^ th mode of the final EMF state. The transition amplitude as written 
above is a functional sum over all distributions {m^}. Differential equations for the functionals that appear in the 
transition amplitude can be found from the finite difference equations of Eqs. (|f 5116(1 . 

F[{ms}} = -^m rt F[K}]+i^G p [{m r i fl }] (25) 

q Ip 

G p [{m^}} = -i(uj + m q uj q )G p [{m i }} + i\ p F[{m^ + 5^}} 

(26) 

*'[{m f }] = -i(uo + ^m^ q )* f [{me}]+i^\ p *f[{m ( :}] 

q p 

+ i$^m,Aj¥f[{m 4 -%}] (27) 
i P 

= -i(H + i;m,w,)^[{mf}]-i5]A,^[{m{}] 



iX p ^ f [{m s + 6 ip }] (28) 
^ [{"»«}] = -*(^ + X)ro9W g )$£[{ro c }]+*Xptf / [{m c }] 

+ ^m g A^[{m £ -5 Sg }] (29) 

+ ^A p $f [{m £ + 5 Cp }] (30) 

Although the transition amplitude of Eq. (|24|l and the differential equations with all Grassmann variables removed 
of Eqs. (|25I30I) can be used from this point onward, it is simpler instead to work with Eq. I|17|l during the trace over 
final EMF states. In the next section we shall combine forward and backward versions of the transition amplitude to 
construct the reduced propagator. 

B. Reduced Propagator 

The evolution of the reduced system with an initial atomic state is given by, 

P (t) = / dfj,(Tn)dfi(r]i) Y[ [dfj,{z ik )dfi(z' ik )} J R (t, Q)iZ(O), (31) 

J k 

in which i?(0) is the combined initial system-bath density operator and Ju(t, 0) is the propagator for the reduced 
system, 



Mt,0) = J dn({z fk })K(t,0)K'{t,0). (32) 
Carrying out the integration with Eq. (|17|l and its barred conjugate one finds, 

J R (t, 0) = exp{fj f m + $'(tW f + Vf9k(t) + 9k(tW f 

ft ft , 33 . 

+ + &(*))(/*(*) + &(*))} 
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C. Initial Thermal State 

For thermal vacuum the initial state in the coherent state representation and in units such that Boltzmann's 
constant is unity (k^ = 1) is, 



R(0) 



J|exp{ 



6 ^ ZifcZ,- 



fc} 



x [poo + ViPlO + ViPOl + ViViPU 



(34) 



Evaluating Eq. (|3 lf> with substitutions from Eq. (|33[) and Eq. <|34|) one may obtain the evolved reduced density 
operator. After expanding completely, the reduced density matrix elements become, for the upper state population, 

Pn(t) = Poo £ ^m^m^ - %}]G;[{m c - %}] e ~^« m ^ 

+pu E ^hh+E^^ik-^}] 

{m £ } V I / 



(35) 



for the lower state population, 



{■m 5 } i 



and for the off-diagonal, 



Pio(i) = pio 



Pm F[{m^}F'[{m ( }\ 



{m k } \ I / 



m q u) q 



(36) 



(37) 



in terms of the definitions of Eqs. (|18I23(I . with {pn, pio, Poi, Poo} being the initial qubit density matrix elements. 



1. Low temperature 

The computation of the reduced density matrix elements involves the calculation of the functionals of Eqs. i|18l 
12311 and the evaluation of the functional summations in Eqs. i|35l37[l . In order to calculate the functionals, a low 
temperature and a weak coupling approximation are applied to Eqs. l|25l3Uf) . Details of the calculation are shown in 
Appendix A. The resulting expressions for the reduced density matrix elements, valid at low temperature (e~^ Wo << 1) 
and weak coupling (A 2 << 1), are 



PiiW = 



1 - T(t) 



Poo 



1 - 



1 - e- r °* 



Pit 



Poo(*) = T(£)poo 



1 - e 



-r„t 



TWpi 



with the definition 



T(t) 



1 - e-P u 



(38) 

(39) 
(40) 

(41) 



and r o = 2A "° being the zero temperature spontaneous emission rate. These reduced density matrix elements are 
illustrated in Fig. (1). 



FIG. 1: These plots illustrate the non-Markovian reduced qubit matrix elements from Eqs. 138I41H for the case of an initial 
a x eigenstate, (|0) + |l})/\/2 at low temperature (e _/3 "° = 0.05), versus a dimensionless time in units of T^ 1 where r o is 
the zero temperature spontaneous emission rate. The diagonal matrix elements thermalize to the low temperature values of 
Eqs. 1)42143^ . and the off-diagonal matrix elements decohere non-exponentially. 

In the long time limit the populations tend to the following thermal values valid at low temperature, 

Pll (t^oo) = e-P»° (42) 
poo(*->oo) = l-e-^° (43) 

and the off-diagonal coherence decays completely 

pio(i^oo) = 0. (44) 

2. Zero temperature limit 

At zero temperature (3 = oo and Eqs. (|38I4()(I become, 

Pii(t) - Pne~ rot (45) 
PooW - Poo + Pn (l - e- r °^ (46) 
Pio(*) = Pioe- r ° t/2 - iw ° 4 (47) 
which is the expected result from Ref. [Tl| . 

3. The Markov approximation limit 

For reference purposes we also include the results in the Markov approximation, which are valid in the regime of 
high temperature which are, 

Pn(t) = Pii(0)e- r ° coth(/W2)t + - - e _^ o (1 - e - r ° coth (^°/ 2 )*) (48) 

poo(i) = l-pu(t) (49) 
PioW = pi (0)e^ t - E ^ ooth < /Sw "/ a >. (50) 



Their asymptotic values are, 



Pll (t - oo) = T^^T « e ~^° + O((o-^) 2 ) (51) 
Pm {t -> oo) = 1 + l_ 0Uo « 1 - e-^- + 0((e-^) 2 ) (52) 
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FIG. 2: This plot shows the quantity , (the ratio of the non-Markovian decoherence rate defined by Fdec (t) :— — , over 

the zero temperature decoherence rate, r o /2) versus the dimensionless time T t in the low temperature regime, e -13 ^" = 0.05. 
The dotted line is the value of the Markovian prediction at finite temperature, coth(/3u;o/2). Initially the non-Markovian 
prediction matches the Markovian result (at the dotted line). As the qubit and EMF become correlated the reduced dynamics 
deviates from the Markovian prediction and the decoherence rate asymptotes to the zero temperature value (along the dashed- 
dotted line). 

and the off-diagonal coherence decays completely 

p 10 (i^oo) = 0. (53) 

The thermal populations in the non-Markovian low temperature approximation match the Markovian thermal popu- 
lations up to 0((e^°) 2 ). 

IV. DISCUSSION 
A. Decoherence 

The decoherence rate is found by computing the off-diagonal elements of the reduced density matrix (e.g. pio(i)). 
The inclusion of bath as well as system dynamics causes the fall off of the off-diagonal matrix elements to become 
slightly sub-exponential. From previous work |nj we know that at zero temperature the decoherence rate is To/2 = 
X 2 uj /tt. Markovian approaches (e.g. 01) predict a decoherence rate of T coth(/?cJ /2)/2, valid at high temperatures. 
The present calculation shows that the decoherence rate, Td ec {t) ■— — N| , actually changes as the total system 
evolves. As shown in Fig. (2) the decoherence rate at t — 0, when the bath is by assumption in a thermal state 
uncorrelated with the qubit, agrees with the prediction of Markovian approaches. As the system and bath evolve 
together the decoherence rate falls back down to the zero temperature value. The interpretation of this is: initially 
the two cases have the same decoherence rate because by arrangement the combined system is a product state of qubit 
and thermal bath, which is the state assumed in Markovian approaches (there is no prior correlation). As the system 
and bath interact, the correlations that arise alter the reduced system dynamics and the combined state evolves away 
from that initial factorizable state. The overall effect is that the the qubit decoheres more slowly in non-Markovian 
dynamics than in Markovian dynamics. 

B. Relaxation 

The relaxation time scale is measured by the value of pn(t), assuming that pn(0) — 1. Similar to the case 
of decoherence, because the initial state of the combined system-bath is taken to be a product state of qubit and 
thermal bath, as it is in Markovian approaches, the dynamics for the populations initially agree in non-Markovian 
and Markovian dynamics (see the inset of Fig. (3)). Then as the system and bath interact, the non-Markovian result 
for the dynamics of the reduced system, which takes into consideration the dynamics of both the bath and the qubit, 
deviates from the Markovian prediction, as shown in Fig. (3). However, the long time behavior of our prediction 
matches the thermalization prediction of the Markovian prediction up to 0((e _ ^ w °) 2 ). Most importantly, Fig. (3) 
shows that the upper state population relaxes more slowly in non-Markovian dynamics than in Markovian dynamics. 
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FIG. 3: This plot shows the difference in the upper state populations, pn(i) — pri aT ° V (£), between the non-Markovian 

prediction, pn(t), and the Markovian prediction, Pn ar ^ ov (t), given that the qubit is initially in the upper state (i.e. pn(0) = 1). 
The difference is plotted versus dimensionless time Y t and in the low temperature regime, e _/3 "° = 0.05. Inspection of the plot 
shows that in non-Markovian dynamics the upper state decays more slowly than in Markovian dynamics. At long times the 
difference in the populations is zero up to 0((e - ' 3 " ) 2 ) (see Eq. I|51|l '). The inset shows that the non-Markovian and Markovian 
predictions agree initially. 
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FIG. 4: This plot shows the ratio of the non-Markovian decoherence rate to the non-Markovian relaxation rate as a function 
of dimensionless time, F t, in the low temperature regime, e~^"° = 0.05. The value is approximately constant at 0.5, which 
is also the ratio between the Markovian decoherence and relaxation rates. The is because both cases share the same physical 
determining factor, i.e. that the resonant type of coupling is at work in this model. 



We may define the relaxation rate (for the case that pn(0) — 1) as T re i(t) :— — ^7^3^~r^ ■ The plot of Fig. (4) 

shows the dependence of the ratio Tdec(t)/T re i(t) on time. It demonstrates that the relaxation and decoherence rate 
are of the same order of magnitude. In other words, the rate of quantum phase information escaping from the system 
to the environment is the same as the rate of energy flow. This property is characteristic of the resonant coupling 
between the two-level atom and the EMF, which leads to a different decoherence behavior from quantum Brownian 
motion (QBM) models One way to visualise the distinction is the realization that in QBM the couplings allow 
the interaction of the system with the far- infrared modes of the environment. The system then loses the phase 
information through soft photons which however carry very little energy. Hence in QBM systems, the relaxation time 
is much longer than the decoherence time. However, in resonant systems, such as being studied here and in the 
system interacts primarily with the modes of the environment near the resonance frequency. Consequently, the phase 
information escapes through photons of energy equal to that of the atom and the decoherence rate is essentially the 
same with relaxation rate. We should remark that although the present results only hold for the low-temperature 
limit the near equality of decoherence and relaxation rate is valid even in the high temperature limit as can be seen 
already from the Markov approximation. 
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FIG. 5: This plot shows the quantity f = Tr[p{t)U o (t)p(0)U^(t)] as a function of the dimensionless time, T t, and at low 
temperature (e~ /3 "° = 0.05) for an initial a x eigenstate, (|0) + |l})/\/2), with U (t) being the free evolution operator. Being a 
measure of the persistence of the initial qubit state after interaction with the environment, it can be considered as the fidelity 
of the qubit in its environment. The non-Markovian fidelity is plotted as a solid line and the Markovian fidelity is plotted 
as a dashed line. The inset is the non-Markovian fidelity minus the Markovian fidelity. Inspection of it shows that in the 
non-Markovian dynamics the EMF bath degrades the fidelity of the qubit more slowly than in the Markovian dynamics. 
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FIG. 6: This plot shows the von Neumann entropy, s(t) = — Tr plog e (p), of the reduced qubit density matrix versus the 
dimensionless time, T t, for low temperature (e~' 3 " = 0.05). The von Neumann entropy is a measure of the purity of a 
density operator. Both the non-Markovian and Markovian von Neumann entropies are plotted, as solid and dashed lines, 
respectively. Inspection of the plot shows that as the qubit interacts with the environment it becomes more mixed, but as the 
qubit equilibrates with the environment (after a time on the order of the relaxation time) it becomes less mixed due to the 
low temperature. In the inset the difference of the non-Markovian von Neumann entropy minus the Markovian von Neumann 
entropy is plotted. It shows that during the initial period of mixing non-Markovian dynamics predicts less mixing than does 
Markovian dynamics. Then, during the later equilibration period Markovian dynamics predicts a less mixed state. Finally, 
both dynamics reach a low temperature equilibrium state which is less mixed than at intermediate times. As in Fig. (3), the 
non-Markovian dynamics asymptotes to thermal equilibrium more slowly than the Markovian dynamics. 



C. Entanglement 

There exists no computable measure of entanglement between a qubit and an infinite continuous bath such as the 
electromagnetic field. However, since Markovian predictions explicitly exclude system-bath entanglement, comparison 
of those predictions with the present results can reveal the effects of system-bath entanglement. First is the deco- 
herence rate discussed previously (see Fig. (2)). Its evolution from the thermal to the zero temperature value shows 
that the combined system-bath reaches and holds at some state in which the off-diagonal elements of the qubit are 
no longer affected by the thermal nature of the bath. No product state could give such behavior and cause thermal- 
ization of the populations. Second is the fidelity, f = Tr[p(t)U o (t)p(0)U%(t)}, of the (|0) + \l))/y/2 state shown in 
Fig. (5), with U (t) being the free evolution operator. The fidelity in this case is a measure of the persistence of the 
initial qubit state after interaction with the environment. Inspection of the difference between the non-Markovian 
and Markovian predictions for fidelity (inset of Fig. (5)) shows that non-Markovian dynamics predicts a slower loss 
of fidelity than Markovian dynamics, although with continued interaction both dynamics predict complete loss of 
fidelity. Third is the von Neumann entropy, s(t) = — Tr plog e (p), for the initial qubit state |1), shown in Fig. (6). 
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The von Neumann entropy is a measure of the purity of a density matrix. Inspection of the difference in this case 
(inset of Fig. (6)) shows that the Markovian result initially predicts a greater loss of purity, but after a time on the 
order of the decay timescale, it predicts less loss of purity than the non-Markovian result. The reason for this seeming 
contradiction is that at low temperature, the system in its approach to thermal equilibrium, is driven to a less mixed 
state. Since correlations with the bath slow the drive to thermal equilibrium in non-Markovian dynamics, as shown in 
Fig. (3), this process is slower in the non-Markovian regime. Comparison of these three quantities (decoherence rate, 
fidelity and von Neumann entropy) shows a consistent picture in which non-Markovian dynamics is characterized by 
the preservation of coherence for longer time. 



D. Conclusion 



We have studied a two level atom coupled to an electromagnetic field (EMF) at finite temperature in the multimode 
Jaynes-Cummings model. We have computed the reduced evolution of the two level system and addressed the issues 
of decoherence, relaxation, and entanglement from its interaction with the EMF bath. Our approach makes use of a 
modified influence functional technique, which enables one to compute the reduced system dynamics while including 
the entangled evolution of the bath and qubit degrees of freedom. That is in contrast to Markovian approaches, 
which assume a fixed bath and hence neglect any dynamics in the bath. We adopt a Grassmann coherent state path 
integral representation for the atom degrees of freedom and bosonic coherent state path integral representation for the 
electromagnetic field, and assume a weak coupling (2 nd order Born) approximation under low temperature conditions. 

We find non-exponential decay for both the diagonal and off-diagonal matrix elements of the qubit 's reduced density 
matrix. From the diagonal elements we see the qubit relax to thermal equilibrium with the bath. From the off-diagonal 
elements, we see the decoherence rate beginning at the rate usually predicted for a thermal bath, but evolving to 
the zero temperature decoherence rate as the qubit and bath become entangled. Comparison of the relaxation and 
decoherence rates, shown in Fig. (4), reveals that as in the zero temperature case, both rates are comparable. At short 
times the ratio of the decoherence to the relaxation rate is initially smaller, but only by a small amount. At higher 
temperatures the initial difference between the two does increase, but that regime reaches the limits of validity of the 
present results. We can see why at low temperatures both rates are related to the atomic transition rate, because it 
is the only relevant physical scale present (unlike at finite temperature where the thermal scale is also at work). This, 
in turn, is a consequence of the particular resonant coupling between the two-level atom and the EMF, as explained 
earlier in 11]. 

On the entanglement of a qubit with the EMF (under resonant coupling) we calculated the qubit 's fidelity and 
the von Neumann entropy. The Markovian result predicts higher loss of fidelity and purity as compared with the 
more accurate non-Markovian calculations. Qualitatively, with the inclusion of quantum correlations between the 
qubit and its environment, the non-Markovian processes tend to slow down the drive of the system to equilibrium, 
prolonging the decoherence and better preserving the fidelity and purity of the system. 

Acknowledgements This work is supported in part by ARDA contract MDA90401/C0903, a NSF and a NIST grant 
to the University of Maryland. 



E. Knill, R. Laflamme and W. Zurek, Proc. R. Soc. Lond. A 454 365 (1997); E. Knill, Keynote address given at First 
International Quantum Information and Quantum Control Conference, Fields Institute, University of Toronto, July 19-23, 
2004 ( http: / /www.fields.utoronto.ca/programs /scientific/04-05 / quantumlC/abstracts /knill. pdf ) . 
M. Abdel-Aty, S. Furuichi and S. Nakamura, Quant. Inf. Comp. 2 272 (2002) 
J. Eisert, C. Simon and M. B. Plenio, J. Phys. A: Math. Gen. 35 3911 (2002) 

M. Keyl, D. Schlingemann and R. F. Werner, Quant. Inf. Comp. 3 281 (2003) 

G. K. Brennen, Quant. Inf. Comp. 3 619 (2003); G. K. Brennen and S. S. Bullock, Phys. Rev. A quant-ph/0406064 
J. I. Latorre, E. Rico, and G. Vidal, Quant. Inf. Comp. 4 48 (2004) 
J. P. Barnes and D. J. Warren, Phys. Rev. A. 60 4363 (1999) 
J. Gea-Banacloche, Phys. Rev. A 65 022308 (2002) 
S. J. van Enk and H. J. Kimble, Quant. Inf. Comp. 2 1 (2002) 
A. Silberfarb and I. H. S. Deutsch, Phys. Rev. A 69 042308 (2004). 



C. Anastopoulos and B. L. Hu, Phys. Rev. A 62 033821 (2000) quant-ph/9901078 

S. Shresta and B. L. Hu, Phys. Rev. A 68, 012110 (2003) \gumt-ph/0301180\ 

S. Shresta, B. L. Hu and N. G. Phillips, Phys. Rev. A 68, 062101 (2003) quant-ph/0302004 

H. J. Carmichael Statistical Methods in Quantum Optics 1: Master equations and Fokker-Planck Equations (Springer 
Berlin, 1999). 



12 



[15 
[16 
[17 

[is; 

[19 
[20 
[21 

[22 



[23 
[24 

[25 
[26 
[27 
[28 
[29 
[30 
[31 

[32 
[33 
[34 



W. H. Louisell, Quantum Statistical Properties of Radiation (J Wiley and Sons, 1973). 
D. F. Walls and G. J. Milburn, Quantum Optics (Springer Verlag, Berlin, 1994). 

C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, Atom-Photon Interactions: Basic Processes and Applications (J 
Wiley and Sons, 1992). 

0. Scully and M. Suhail Zubairy, Quantum Optics (Cambridge University Press, Cambridge, 1997). 

C. W. Gardiner and P. Zoller, Quantum Noise (Springer Verlag, Berlin, 2000). 

R. Alicki, M. Horodecki, P. Horodecki, and R. Horodecki, Phy. Rev. A 65 062101 (2002). 

J. I. Kim, M. C. Nemes, A. F. R. de Toledo Piza and H. E. Borges, Phys. Rev. Lett. 77 207 (1996). 

H. -P. Breuer, A. Ma and F. Petruccione, Time-local master equations: influence functional and cumulant expansion 
Quantum Computing and Quantum Bits in Mesoscopic Systems eds. Anthony Leggett, Berardo Ruggiero, Paolo Silvestrini 
(Kluwer Academic Plenum Publishers, 2003) quant-ph/0209153 
M. Esposito and P. Gaspard, Phys. Rev. A 68 066112 (2003). 

D. Tolkunov and V. Privman, Phys. Rev. A 69 062309 (2004); V. Privman, J. Stat. Phys. 110 957 (2003); V. Privman, 
Mod. Phys. Lett. B 16 459 (2002). 

J. Gambetta and H. Wiseman, Phys. Rev. A 66 052105 (2002). 

T. Yu, quant-ph/0402086 

F. Haake, Z. Physik 223 353 (1969). 

N. G. van Kampen, Stochastic Processes in Physics and Chemistry (Elsevier Science B. V., Amsterdam, 2003). 
S. Shresta, Ph. D. dissertation, University of Maryland 2003. 

1. M. Gelfand and A. M. Yaglom, J. Math Phys. 1 48 (1960). 

Y. Ohnuki and T. Kashiwa, Coherent states of Fermi operators and the path integral Coherent states: applications in 
physics and mathematical physics eds. J. Klauder and B. Skagerstam (Singapore: World Scientific 1978) 449-465. 
A. Perelomov, Generalized coherent states and their applications (Berlin: Springer 1986). 
J. Seke and W. N. Herfort, Phys Rev. A 38 833 (1988) 
D. P. DiVincenzo and D. Loss, \conlTm at/0405525\ 



APPENDIX A: CALCULATIONAL DETAILS OF QUBIT IN A THERMAL BATH 

1. Approximated functional solutions 

Eqs. H25I30|I are two sets of coupled differential equations. One set being the pair of equations 

F[{ m ii}} = -i^m q u q F[{Tn i }}+i^mi\iG p [{m i '' S^i}} (Al) 

q Ip 

G p [{m^\] = -i(uj + ^m q uj q )Gp[{m i }} + i\ p F[{m^ + <%,}] (A2) 
q 

and the remaining four equations comprising the other set. The solution method for this pair in the low temperature 
and weak coupling limits will be sketched out in this appendix. The solutions for the other set in the same limits will 
follow a similar sequence. First, given the initial conditions 

F[{m e }](t = 0) = 1 (A3) 
G p [{m $ }](t = 0) = (A4) 

the Laplace transforms of the above equations are 

zF[{m ( }]{z) - 1 = -i^mqUqFlim^W + i^miMGplimt-Szijjiz) (A5) 

q ip 

zG p [{m{\](z) = —i(u + y^ j m q oj q )Gp[{m^}](z) + iXpF[{m^ + S^ p }](z). (A6) 

q 

The second equation can be rearranged into 

GA{rnMz)= iX :f^ P}]{Z) V ^ 
z + i(lj + m q uj q ) 

which can be substituted back into Eq. (|A5|I to give 

F[{m,}]( Z ) = i ^y^^/N-^y;^ (A8) 
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In the expression above the low temperature approximation is applied by setting p — I in the summation of the RHS. 
The justification is that the summation on the RHS will be peaked about lui = lu such that the greatest contribution 
from F[{m^ — 5^i + 5^ p }](z) will be for ui = uj . However, at low temperatures those frequencies will not be populated. 
As a result the vacuum will be annihilated, unless S^ p — S^i, which will cause the major contribution from the p 
summation to be from p = I. The low temperature approximation is thus that the temperature is low enough that the 
modes with frequency oj are unoccupied, i.e. e _ ^ Wo << 1. Applying this approximation, Eq. I|A8J) can be rewritten 
as 



F[{ro e }](z) = iz + i^2m q uj q + Y^ 



TO/A? 



q i 



z + i(uj -Lu t + J2 q m q u q ) 



The zero order pole of F[{m^}](z) is at , 



The reaction term at this point is found equal to 



r m 



(A9) 



-*A, 



with F = x ^° , which shows that the second order shift in the pole includes both a real and an imaginary part. After 
absorbing the imaginary part in a renormalization of the frequency, the second order pole is z = — i £ m„uj„ ' 
with the definitions 



E 



TO, 



(A10) 



The desired functional can be calculated as in inverse Laplace transform of Eq. (| A9|) at the second order pole to give 



'. m q UJ q t 



(All) 



The inverse Laplace transform contains a contribution of a branch cut as well as a pole |33j, |3^|. We ignore the 
contribution of the branch cut, which is negligible at all but very late times such that T a t > 20 and very early times 
such that T t < 1CP 21 (see Eq.(3.20) of Ref. 33]). In all cases, we assume that time is much later than the inverse 
cut-off time. Further comparison between the branch cut and the non-Markovian correction over Markovian dynamics 
shows that the branch cut contribution is smaller by greater than three orders of magnitude for T t > 0.1. The other 
functional in the pair can be calculated by integrating Eq. (|2*B|) 



Gi[{mt-5 & }](t) = i 



A 1 



'-t—iJ2„ m q UJ q t 



r„ 



+ i(uJi - LU ) 



Following similar calculations the rest of the functionals are found to be 

. r (m + i) t _i( 0Jo +Y l m q ui q )t 



* / [{mJ](t)=e" 
^[{mc}](t)= A ' 



Xe -^t-i(u, +E 9 



,-r m 



UJp{u) p — LU 
'-t—i{w„+J2 „ m q u q )t 



UJ p {uJp — L0 o — i-f) 

^g-i^o+Ui-Up+E, m q uj q )t 



/UlyJUJp 



1 1 



L t — i(uj p — L0 o )t 



[(wj - uj ) + i^f] [(w p - lo ) - i- 



(A12) 

(A13) 
(A14) 

(A15) 



(A16) 



2. Computation of density matrix elements 

The solutions of Eqs. I|A11IA16|) can be substituted into Eqs. (|35I37|) to evaluate the reduced density matrix 
elements in the limits of low temperature and weak coupling. The reduced density matrix elements in that form are 
summations over all distributions {m^}. The pw(t) matrix clement will be demonstrated below as a representative 
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calculation. The evaluation of the other summations follow along similar lines. From Eq. (|37|l . the off-diagonal density 
matrix element is 



{m ? } V ; / 



First, from Eqs. (1 A 111 A 161) the functional in parentheses can be determined to be 



(A17) 



(A18) 



so that the off-diagonal matrix clement becomes 



Pw 



M sr J r (2m + i) 
(*) = pw E exp \ 2 



(A19) 



Denoting by primes those terms for which lu^ — uj and double primes those for which lu^ =^ lu , the summand can be 
rewritten with the substitution ra Q = ro£, 



Pw(t) = Pw 



{™ 5 } « £ 



(A20) 



The summation over distributions can be more clearly written as 



E 

{m £ } 



IIE 

^ m 5 =0 



nE 



which allows us to bring Eq. I|A20|) into the form 



Pw(t) = Pw e 2 



1 - 



1 _ g-(r o t+/3w ) 



nE 

5 m ? =0 



(A21) 



(A22) 



The factor at the end is removed by normalization of the reduced matrix element by its value if r o = 0. The final 
result for the off-diagonal matrix element is 



Pw(t) = Pw e" 



1 - e~ Pu 



1 _ e -(r t+i3uj ) 



(A23) 



with r o being the zero temperature spontaneous emission rate. The other reduced density matrix elements are given 
in the text. 



